home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / dpotrf.z / dpotrf
Text File  |  1996-03-14  |  3KB  |  133 lines
  1.  
  2.  
  3.  
  4. DDDDPPPPOOOOTTTTRRRRFFFF((((3333FFFF))))                                                          DDDDPPPPOOOOTTTTRRRRFFFF((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      DPOTRF - compute the Cholesky factorization of a real symmetric positive
  10.      definite matrix A
  11.  
  12. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  13.      SUBROUTINE DPOTRF( UPLO, N, A, LDA, INFO )
  14.  
  15.          CHARACTER      UPLO
  16.  
  17.          INTEGER        INFO, LDA, N
  18.  
  19.          DOUBLE         PRECISION A( LDA, * )
  20.  
  21. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  22.      DPOTRF computes the Cholesky factorization of a real symmetric positive
  23.      definite matrix A.
  24.  
  25.      The factorization has the form
  26.         A = U**T * U,  if UPLO = 'U', or
  27.         A = L  * L**T,  if UPLO = 'L',
  28.      where U is an upper triangular matrix and L is lower triangular.
  29.  
  30.      This is the block version of the algorithm, calling Level 3 BLAS.
  31.  
  32.  
  33. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  34.      UPLO    (input) CHARACTER*1
  35.              = 'U':  Upper triangle of A is stored;
  36.              = 'L':  Lower triangle of A is stored.
  37.  
  38.      N       (input) INTEGER
  39.              The order of the matrix A.  N >= 0.
  40.  
  41.      A       (input/output) DOUBLE PRECISION array, dimension (LDA,N)
  42.              On entry, the symmetric matrix A.  If UPLO = 'U', the leading N-
  43.              by-N upper triangular part of A contains the upper triangular
  44.              part of the matrix A, and the strictly lower triangular part of A
  45.              is not referenced.  If UPLO = 'L', the leading N-by-N lower
  46.              triangular part of A contains the lower triangular part of the
  47.              matrix A, and the strictly upper triangular part of A is not
  48.              referenced.
  49.  
  50.              On exit, if INFO = 0, the factor U or L from the Cholesky
  51.              factorization A = U**T*U or A = L*L**T.
  52.  
  53.      LDA     (input) INTEGER
  54.              The leading dimension of the array A.  LDA >= max(1,N).
  55.  
  56.      INFO    (output) INTEGER
  57.              = 0:  successful exit
  58.              < 0:  if INFO = -i, the i-th argument had an illegal value
  59.              > 0:  if INFO = i, the leading minor of order i is not positive
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. DDDDPPPPOOOOTTTTRRRRFFFF((((3333FFFF))))                                                          DDDDPPPPOOOOTTTTRRRRFFFF((((3333FFFF))))
  71.  
  72.  
  73.  
  74.              definite, and the factorization could not be completed.
  75.  
  76.  
  77.  
  78.  
  79.  
  80.  
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.